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Abstract 

We present an ansatz for the quark and lepton mass matrices, 
derivable from SO(IO) type GUTs, which accommodates a heavy (> 
92GeV) top quark and permits large mixings in the <-> Vr sector (as 
suggested by the recent Kamiokande and 1MB data on the atmospheric 
neutrinos). The well known asymptotic relations = m,-, = ^fn^ 
and rridms = niem^ all hold to a good approximation. Depending on 
mixing which can even be maximal, the mixing angle relevant 
for solar neutrino oscillation lies in the range 7.8 x 10~^ ;S sin^ 20e^ ~ 
2.1 X 10^^. For the ^^Ga experiment the event rate, normalized against 
the standard solar model prediction of 132 SNU, is estimated to be 
between 80 and 20 SNU. 
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There is currently a surge of experimental activity around the excit- 
ing possibility that one or more species of the neutrinos may possess a 
tiny mass. The most recent analysis*^^^ of the Kamiokande^^^ and the 1MB 
data*^'^) on atmospheric neutrino interactions suggest a discrepancy with the- 
oretical expectations*^'^'' which can be nicely explained in terms of neutrino 
oscillations^-'^'^) . Assuming a two flavor (z/^ ^ u^-) oscillation, the relevant 



mixing angle satisfies the constraint sin 29^r ~ 0.42. The mass difference 
squared lies in the range (10~^ — 10^^')eV^. 

Perhaps the most surprising aspect of this observation, if true, is the rel- 
atively large vahic suggested for the — mixing angle 9^r{^ 20°). Within 
the framework of grand unified theories with the minimal higgs structure, 
for instance, one typically finds that^^'''^ O^r ~| Vcb \— 0.05, where Vcb is the 
(c6) element of the well known KM matrix of the quark sector. By making 
some minor modifications in the higgs structure, one could perhaps make O^r 
three times larger than | |, but this still is much too small compared to 
what the atmospheric neutrino data suggests. 

The purpose of this paper is to address this and some important re- 
lated issues within a more general framework, which is inspired by GUTs as 

well as some earlier work of Pritzsch^^^ and others. The idea is to write 
down an ansatz for the fcrmion mass matrices which i) is predictive in 
the quark sector and admits a heavy (> Si2GeV) top quark, ii) preserves 
to a good approximation some well known asymptotic relations including 

rri}, ~ rriT-jms — ^m^^mdms — mem^ and 6c — {md/ms)^, and iii) allows for 
large mixings in the z/^ ^ Ur sector. It turns out that in the quark sector 
the new ansatz we propose has essentially the same predictive capacity as 
the original one of ref. (8). A particularly important feature of the new 
scenario is the restriction on the mixing angle relevant for the solar neutrino 
problem. One finds a significant deficit relative to the standard solar model 
prediction, which will soon be tested in the ongoing SAGE^^^ and GALLEX 
experiments. 

The ansatz that we will consider can be motivated within a grand unified 
framework (e.g., SO(10))(i°) with a non-minimal higgs structure. It works 
both with supersymmetric (SUSY) as well as with non-SUSY GUTs. In 
the latter case, in order for the proton lifetime to be compatible with the 
experimental lower bound, an intermediate step (e.g., SU{A)c x SU{2)l x 
SU{2)r)''^^^ would be needed. Consider then the following ansatz for the 
quark and lepton mass matrices: 




Mh = a De'" B I M„ 
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Several comments are in order: 

1. Wc have written down the matrices of Eq. (1) after a suitable redefini- 
tion of the phases of the fermion fields. The non-zero elements of the 
matrices, before this redefinition, are allowed to have arbitrary phases. 
All but one phase a can be removed from the Dirac mass matrices. 
The charged current interaction in this basis is not proportional to the 
identity matrix, but has the generation structure given in terms of two 
phase parameters a and r by 



The entries A, A', 5, B\... can be chosen to be real and positive without 
loss of generahty. 

2. Within an SO(IO) type framework, the entries A,A\C^C' arise from 
the higgs 10 plets, while B, B' and D arise from the 126 plets. The 
reason for the form of the D terms is to retain, to a good approximation, 

the well-known asymptotic relation ~ The B,B' terms 

are motivated by the desire to accommodate both large mixings in the 
2-3 lepton sector and a moderately heavy top quark (mj < IQOGeV). 
When the parameter B is set to zero, our ansatz reduces to that of Ref. 
12, a detailed analysis of which is presented in Ref. 13; 

3. With C 3> A, 5, D, one would recover nih = rrir ^^^^ to a good approx- 
imation. However, a small violation induced through mixing of this 
asymptotic equality plays an important role in our analysis; 

4. The form of the 2-3 sector of the heavy Majorana mass matrix is di- 
rectly related to the corresponding sectors of the quark and the charged 
lepton matrices. In other words, the same 126 plet of higgs contributes 
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to these three sectors. The (22) elements D and M2 of Eq. (1) arise 
from the same Yukawa couphng to a 126 of Higgs. Similarly, B, B' 
and M3 result from a common Yukawa coupling to another Higgs 126 
plet. (Note that the 126 generating (22) elements D and M2 should 
be distinct from the 126 that generates (23) elemnt B, B' and M3, if 
it were the same, a non-zero (33) element will also be induced. Conse- 
quently, there is no relation of the type M2/M3 = D/B.) We need an 
additional (independent) contribution in the Majorana sector (the 1-1 
entry) to give a large mass to the heavy VeR, 

5. Since the charged fermion sector is described by 11 fundamental pa- 
rameters including tan/5, the ratio of the two Higgs vacuum expectation 
values, we will have 3 predictions. Additionally, since some of these pa- 
rameters are phase angles whose moduh are constrained to be less than 
unity, we also have approximate relations among the physical observ- 
ables (see below). 

It should be emphasized that the above form of the mass matrices is pre- 
scribed at the GUT scale Mu- 

The real symmetric matrix is diagonalized by an orthogonal transfor- 
mation, OuMuOu — M*"f (tilde denotes transpose), where 



Here m^, —rric and denote the eigenvalues of and we have used the 
relations A' ~ ^/rn^JfTc, B' ~ ^/rnjfrt and C c:^ mt. 

Since and are complex symmetric matrices, their diagonalization 
is achieved via bi-unitary transformations of the type UdMdUd = M*"^, and 

iltMfUi = Mf^""^ . Denoting the eigenvalues as (m^, — m^, m^) and (me, — m^u, m^) 
respectively, the appropriate unitary matrices are: 
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Here we have defined 
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assumed C ^ B,D ^ A, and used the relations 
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The phases 0, are determined in terms of B, C, D and a 
relations 
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Prom the relation > 0, one finds that 



1 m 



^ + > 36^ 
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Employing | cos a |< 1, one can show that 



We therefore see that in the hmit of a strict equahty mt, = rrir, we have 
I I = 3 I I . In practice, e will be a small positive quantity which would 
lead to small violations of these relations. 

In the quark sector the Kobayashi-Maskawa matrix is given by 

Vkm^Ou^ jt/* (11) 

which leads to the following (asymptotic!) matrix elements: 




In addition, the reparameterization invariant CP violating quantity J is given 

by(15) 



J = Im{V^sVci>V:,V:,) 

- [s/^^/^"N^('^ - r) + sin (0 + r - 2a)] 




(13) 



sin (0 — a) 



^sin(0-^)]- 



It is clear from Eqs. (12-13) that once the charged fermion masses are 
specified, all the mixing angles as well as J are determined in terms of the 
phases a and r. This situation is analogous to the Fritzsch ansatz, except 
that our modified version can accommodate a heavy top quark (due to the 
difference in the expression for K&D- 

As was pointed out earlier, the above relations hold at the GUT scale and 
so comparison with data requires that we consider their evolution with mo- 
menta. It was noted quite some time ago^-*^^^ that as far as the mixing angles 
are concerned, significant evolution will occur for the entries involving the 
third family provided that the top quark is sufficiently heavy (> lOOGeV). 

For definiteness, we focus on the supersymmetric case. We further sim- 
plify our analysis by assuming that the sparticles and the second higgs dou- 
blet are degenerate at 300 GeV. Employing the one loop renormalization 
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group equations, and fixing the top quark mass at ISQGeV, we first deter- 
mine the unification scale Mu. As low energy inputs we use the precisely 

known values of ai and eta at M^^^^: ai{Mz) = 0.01013, a2(M2) = 0.03322. 
We find the unification scale to be Mu = 9.8 x lO^^GeV^ and the common 
gauge coupling at Mu to be au{Mu) = ^-g. Running backwards, we deter- 
mine asi^Mz) to be 0.105, which is in agreement with measurements. 

The evolution equations for the elements of the KM matrix can be found 
in ref. (16). It turns out that for tan/? < 7(tan/3 = V2/V1, the well-known 
ratio of the two vevs present in the minimal SUSY extension of the standard 
model), the coupled equations can be solved semi-analytically. Here we will 
content ourselves by presenting results (Table 1) showing the variation of 
the mixing angles with momentum, for varying values of rrit and tan /3. For 
rrit not too large (< ISOGeV), the variation in | Vcb |, | Vub |, | Vtd \,\ Vts \ 
is < 5%. However, if the top is close to its maximal allowed value of about 
190GeV, the variations in some of the quantities can exceed 20%. In Table 
2 we display the dependence of the quark and charged lepton mass ratios 
m(M[/)/m(/i), on rrit and tan/3. For the d and s quarks, as well as for the 
(e, fi, t) leptons these ratios are (essentially) independent of rrit and tan /3 
and are in the range (0.177-0.179) and (0.656-0.689) respectively. 

From Tables 1-2, we can infer the predictions of our ansatz for the quark 
masses and mixing angles at the weak scale. From the inequality of Eq. (10) 
wc see that the parameter e is bounded by e < 0.10. There is also a lower 
bound on e, e > 0.035^ if we assume that the top mass is around 130 GeV 
(see Eq. (12) for \Vcb\)- This lower bound, however, goes away if the top 
quark is much heavier, near its fixed point value. ""^^'^^ If e is set to zero, one 
recovers the asymptotic relation*^^^'^^'^*'^^'' 

I Feb I- irric/mt)^. (14) 

This relation (Eq. (14)) is disfavored phenomenologically (although not ex- 
cluded), both in the non-sup ersymmctric as well as in the supersymmetric 
case. It requires a top quark in the mass range 180 — 220GeV(175 — 190GeV) 
in the non-SUSY (SUSY) case. Moreover | Vcb \ must be > 0.052^^ to be com- 
pared with the experimentally favored range 0.043 ± 0.009. 

From the constraint 0.035 < e < 0.10 one can deduce the allowed range 
for the 6-quark mass. It follows from Table 2 that mb{mb) hes in the range 
(4.3 — 5.3)GeV (for as{Mz) = 0.105), the larger value corresponding to 
smaller e. This is in accord with the value quoted in Ref. 20, mbirrib) = 
4.25 ±0.1 GeV. (We should point out that nib can be brought further down, 
as in Ref. 13, even if e = 0, provided the top quark mass is around 180 
GeV. Here, however, wc arc interested in a moderately heavy top {rrit ~ 
130 GeV).) Since the asymptotic relation rrib — holds to within 10%, it 
follows from Eq. (10) that | |= | | \ should hold to within 5%. This 
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implies ms{lGeV) = 130 — lAOMeV, in relatively good agreement "^^"^^ with 
observations. From the asymptotic relation mdiris = mem^^^-*, we also have 

a prediction for the d quark mass: = (5.7 — 6)MeV. 

Since the variation of the Cabibbo angle with momentum is negligible, we 
essentially have the Fritzsch prediction for | Vus \ ■ On the other hand, | V^b \ 
can increase by as much as 14% (even for rut < IBOGeV) so that at the weak 
scale, I Vub 1= 0.002 — 0.0037. The CP parameter J naturally comes out in 
the range 10~^ — 10~^ consistent with observations in the K meson system. 

To summarize, in the charged fermion sector, our ansatz is quite predic- 
tive. It reproduces the well-known asymptotic relations m^, = m,-, = 
and mams = rriem^ to a good approximation. The KM angles in the quark 
sector are determined in terms of the fermion masses and two arbitrary phase 
angles. Specifically, there are three predictions in the charged fermion sector. 
One of them is the mass relation m^^mg = mgm^, the other two are mixing 
angle relations of Eq. (12) obtained by eliminating the phases a and r. The 
predictive capacity of our ansatz is as good as Pritzsch's for quark mixing 
angles, but the top quark can be heavy in our case. 

Turning now to the neutral lepton sector, the 3x3 light neutrino mass 
matrix is obtained from the see-saw formula^^^^ 

= M^^^™'=(Mf ''^■°^'^"'^)-iM^^™^ (15) 

where 



orana\~l _ 



\ 
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(16) 



[Note that we are making the assumption that the direct contributions to 
the left handed Majorana neutrino masses are negligible. This is justified 
in a SUSY SO(IO) framework, since the tree-level vacuum expectation value 
of the field which supplies such Majorana masses can consistently be set 
to zero.] Here M is an overall superheavy scale, and ri,r2 denote ratios 
of superheavy masses. The parameter r2 is given by r2 = M3/M2. It is 
reasonable to assume r2 to be of order 1, since M3/M2 ~ B/D which is of 
order 1, (see eq. (7)) times ratio of vev's. We will also discuss the case when 
r2 is not of order 1. The matrix M^y takes the form 
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-3 



1-6, 



'r2e 



(17) 

where P is the diagonal phase matrix P = diag[e'^'^^, 1,1]. Observe that M^, 
has five parameters, viz., ri,r2,0i,02 and M. Note however, that the phase 
(pi appears with a very small coefficient making it an irrelevant parameter. 
The four eff^ective parameters describe six observables: three neutrino masses 
and the three neutrino mixing angles. This leads to two predictions. We shall 
sec below that these two are predictions for mixing angles in terms of the 
neutrino mass ratios. 

It is convenient to diagonalize in two stages. Let M'^ denote the part 
of Miy which excludes the two diagonal phase matrices. We can reduce Ml 
to an effective 2x2 matrix (with zeroes along the first row and column) 

through the transformation UiM^Ui, where t/i is a unitary matrix given by 




(18) 



Note that we are justified in neglecting the term ^^^le**^^ as long as ri is 
of order unity. 

The 2-3 sector of the transformed matrix coincides with that of M^. It 
can be diagonalized by a unitary matrix U2 of the type 
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ce'^p+^^ se'^P-^^ 
—se~^^P~^^ rp-^(p+^) 



ce 



(19) 



where c = cos^ and s = sin^. In particular, the expressions for 6 and p are 
as follows: 
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g /m7 I 3 /^(1 _ 3 g-j<^2) + (1 _ 3,/^r2e+^'^2Vi _ 5 g-^02^ 
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1 - 6^^r2e^<^2 
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Y nit -r^ 1 \l nit nit 



1 - 9,/^ COS 02 + 9^ + 18^r| - 27^ cos 02 



(20) 



(21) 



The light neutrino mass eigenvalues (mass eigenstates are denoted 1^1,2,3) 
turn out to be: 



777.; 



SIM 



m - 81 
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1 - 6,/^r2e*'^2 12^2 I -3,/^+9^r2e*'^2 p^/gmcV 



m 



1^3 



M 



1 - 6,/^r2e^'^2 12 ^2 I -3,/^ + 9^r2e^'^2 12 ^ / gm^ V 



(22) 



We are now in a position to write down the elements of the lepton KM 

rlepto 
KM 



matrix Vt'fJ"". Let = U1U2. Then 



-lAepton rrf 



V 



u. 



(23) 



and one obtains for the off-diagonal elements the expressions: 
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I^3mI ~ 



sin 6 



3ecos^ + sin 66'^^-"+^"') 
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Several comments are in order: 

1. Unlike the quark sector, the mixing angles in the lepton sector do not 
run between Mu and the weak scale since the right-handed neutrinos 
are superheavy. 

2. The mass ratio mj^^/rrii,^ is enhanced by almost two orders of magni- 
tude relative to what one naively expects from the see-saw mechanism. 
This is again related to the factor 3 in the ansatz for the fermion mass 
matrices which, in turn, was motivated by the desire to include large 
mixings. The enhancement will be important when we discuss atmo- 
spheric and solar neutrino oscillations. 

3. The mixing angle relevant for the MSW explanation of the solar neu- 
trino puzzle is given by 



K. 



cos6' 



'nip 



m 



(25) 



Note that as ^ — the expression for | V^^e \ essentially coincides with 
the one derived earlier in ref. (6). If, the recent Kamiokande results on 
the atmospheric neutrino survive the test of time, then 6 is expected 
to lie in the range 18° ^ 6 < 45° (see below). Taking all this into 



account, we estimate that 7.8 x 10"^ < sin^ 2^e^ < 2.1 x 10"^. The 

corresponding | Am^ | is in the range (5 x 10~^ — 2 x 10^^)eK^. For the 

Gallium experiments currently under way the estimated event rate^^^^ , 
normalized relative to the standard solar model value of 132 SNU, lies 
between 80 and 20 SNU. 
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4. A combined fit to the solar neutrino and the atmospheric neutrino data 
requires rriu^ = (1.4 - 2.2) x 10-^eV,m^^ = (0.03 - 0.6)eV. The ratio 
^u2/^uz should then lie in the range (2.3 x 10~^ — 7.3 x 10"^). Prom 
Eq. (22), this constraint implies a lower limit on the parameter r2, 
r2 > 0.8. If we restrict to r2 = (0.8 — 3) so that there is no hierarchy in 
the superheavy masses, the 2-3 mixing angle 9 in the neutrino sector is 
in the range 9 = (18° — 45°). This certainly is in the right parameter 
range for atmospheric neutrino oscillations. For large values of r2, the 



angle 9 approaches 3/2 Jnic/mt ~ 7° (see eq. (20)). Even in this 



case, the atmospheric neutrino deficit can be accommodated, if e is 
near its upper limit of 0.1 (see V^^t of eq. (24)). Note that 9 close 
to 45° gives the largest count (~ 80 SNU) for the SAGE/GALLEX 
experiment. The effective <-> Vj- mixing angle for this case lies in the 

range 0.69 < sin22^^^ < 1. 

In conclusion, the scheme presented here was guided by the desire to 
retain to the extent possible the simplicity of the original Fritzsch ansatz and 
come up with some predictions in the neutrino sector. The SAGE/GALLEX 
experiment should soon be able to tell if we are on the right track! 



Note added: 

After the submission of this paper for publication, GALLEX collaboration 
has announced their result: 83±19±8 SNU. This is in remarkable agreement 
with our prediction (see eq. (25) and the subsequent discussions). 
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Table 1: The evolution factor 



for the quark mixing angles 



cb 



Vub I, 

The 



Vijjmt) 
Vij (mu) 

Vtd \ and I Vts \ (as functions of mt (in GeV) and tan/3), 
corresponding factor for the CP paramaeter J is given by the square of 
these numbers. The running of all the other elements of V are negligible. 
The dashes indicate that the top quark Yukawa coupling has become non- 
perturbative before the GUT scale. 



mt 

itan/? 


100 


130 


160 


190 


1 


1.051 


1.139 






2 


1.027 


1.055 


1.122 




3 


1.023 


1.046 


1.093 


1.359 


4 


1.022 


1.043 


1.085 


1.244 


5 


1.021 


1.041 


1.081 


1.217 


6 


1.021 


1.041 


1.080 


1.205 


7 


1.021 


1.040 


1.079 


1.198 
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Table 2: The evolution factors 



for the quark masses as function 



mi{Mij) 

of rrit and tan /3. For u quark, n = IGeV is chosen and for c and b quarks, 
H = 1.27GeV and 4:.25GeV respectively. The corresponding factors for {d, s) 
and {e,fi,T) are independent of rrit and tan/? and are equal to 0.178 and 
0.673 respectively. 



mt 




1 nn 

lUU 


1 Qn 


1 fin 

iDU 


1 on 
iyu 


i tan/3 
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U.ZUi 


U.ZO ( 






J- 




U.ZiZ 










h 


n 949 


963 








u 


U.io/ 


n one; 


u.Z4y 




2 


c 


u. iyo 


n 01 R 

U.ZiO 


n OR'? 






b 


236 


0.244 


260 






u 


0.185 


0.199 


0.230 


0.448 


3 


c 


0.196 


0.211 


0.243 


0.474 




b 


0.236 


0.242 


0.254 


0.317 




u 


0.184 


0.198 


0.225 


0.343 


4 


c 


0.195 


0.209 


0.238 


0.363 




b 


0.235 


0.241 


0.252 


0.290 




u 


0.184 


0.197 


0.223 


0.322 


5 


c 


0.194 


0.208 


0.236 


0.340 




b 


0.235 


0.241 


0.251 


0.284 




u 


0.184 


0.196 


0.222 


0.312 


6 


c 


0.194 


0.208 


0.234 


0.330 




b 


0.235 


0.241 


0.251 


0.281 




u 


0.184 


0.196 


0.221 


0.307 


7 


c 


0.194 


0.207 


0.234 


0.325 




b 


0.235 


0.240 


0.250 


0.280 
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